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ABSTRACT 



> 

^ Aims. We present full evolutionary calculations appropriate to post-AGB PG1159 stars for a wide range of stellar masses. 

, Methods. We take into account the complete evolutionary stages of PG1159 progenitors starting from the Zero Age Main 
Sequence. We consider the two kinds of Born Again Scenarios, the very late thermal pulse (VLTP) and the late thermal pulse 
(LTP), that give rise to hydrogen-deficient compositions. The location of our PG1159 tracks in the effective temperature - 
fl") . gravity diagram and their comparison with previous calculations as well as the resulting surface compositions are discussed at 
some length. 

Results. Our results reinforce the idea that the different abundances of 14 N observed at the surface of those PG1159 stars with 
undetected hydrogen is an indication that the progenitors of these stars would have evolved through a VLTP episode, where 
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O ' most of the hydrogen content of the remnant is burnt, or LTP, where hydrogen is not burnt but instead diluted to very low 
surface abundances. We derive new values for spectroscopical masses based on these new models. We discuss the correlation 
■ between the presence of planetary nebulae and the 14 N abundance as another indicator that 14 N-rich objects should come from 
*p. a VLTP episode while 14 N-deficient ones should be the result of a LTP. Finally, we discuss an evolutionary scenario that could 
^> , explain the existence of PG1159 stars with unusually high helium abundances and a possible evolutionary connection between 
these stars and the low mass O(He) stars. 

Key words, stars: evolution — stars: abundances — stars: AGB and post-AGB — stars: interiors — 



1. Introduction 

It is accepted that the hot, luminous PG1159 stars consti- 
tute the immediate descendants of post-asymptotic giant 
branch (post-AGB) evolution. They provide an evolution- 
ary link between these stars and most of the hydrogen (H)- 
deficient white dwarf (WD) stars. Specifically, PG1159 
are thought to be formed via a born again scenario, that 
is, a very late thermal pulse (VLTP) experienced by a 
hot WD during its early cooling phase (Fujimoto 1977, 
Schonberner 1979 and Iben et al. 1983 for earlier refer- 
ences) or a late thermal pulse (LTP) that occurs dur- 
ing the post-AGB evolution when H burning is still ac- 
tive (see Blocker 2001 for references). During the VLTP, 
an outward-growing convection zone powered by helium 
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burning shell develops and reaches the H-rich envelope of 
the star, with the consequence that most of the hydro- 
gen content is burnt. The star is then forced to evolve 
rapidly back to the AGB and finally into the central star 
of a planetary nebula at high effective temperatures (T e g) 
but now as a H-deficient, quiescent helium- burning ob- 
ject. LTP also leads to a hydrogen-deficient composition 
but as a result of a dilution episode. Consequently, while 
post- VLTP models are expected to become H-depleted 
(DB) white dwarfs (Althaus et al. 2005a), it has been ar- 
gued that post-LTP models may finally (after the PG1159 
stage) evolve into H-rich (DA) WDs with thin H envelopes 
(Althaus et al. 2005b). 

Currently, 37 stars are members of the PG1159 family, 
which span a wide domain in the log g-T c ff diagram: 5.5 <^ 
log g < 8 and 75000 K £ T cff ^ 200000 K (see Werner & 
Herwig 2006 for a review). Roughly half of them are still 
embedded in a planetary nebula. A striking feature char- 
acterizing PG1159 stars is their peculiar surface chemical 
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composition. Indeed, spectroscopic analyses have revealed 
that most of PG1159 stars exhibit H-deficient and helium- 
, carbon- and oxygen-rich surface abundances. Surface 
mass abundance of about 0.33 He, 0.5 C and 0.17 O are 
reported as typical, though notable variations are found 
from star to star (Dreizler & Heber 1998; Werner 2001). In 
fact, almost every star has its individual He/C/O mixture. 
According to Werner (2001) helium abundances range be- 
tween 33-85 % and carbon and oxygen abundances span 
the ranges 13-58 and 2-24 %, respectively. The rich variety 
of surface patterns observed in PG1159 stars poses indeed 
a real challenge to the theory of post-AGB evolution. In 
particular, the appreciable abundance of oxygen in the at- 
mospheres of these stars has been successfully explained 
by Herwig et al. (1999) on the basis of evolutionary calcu- 
lations of the born-again scenario that incorporate convec- 
tive overshoot. In addition, hydrogen has been detected in 
some PG1159s (called hybrid-PG1159 stars). These H-rich 
PG1159s are currently interpreted in terms of a final he- 
lium shell flash experienced by the progenitor stars imme- 
diately before their departure from the AGB (AGB final 
thermal pulse or AFTP scenario, Herwig 2001). Interest 
in PG1159 stars is also motivated by the fact that about 
ten of them exhibit multiperiodic luminosity variations 
caused by non-radial <?-mode pulsations, a fact that has 
attracted the attention of researchers (see Quirion et al. 
2004; Gautschy et al. 2005; Corsico & Althaus 2006 for 
recent works). 

Needless to say, many relevant aspects about PG1159 
stars, for instance mass determinations, precise asteroseis- 
mological fittings, pulsation stability analysis and ages 
call for the need of extensive full evolutionary models 
of PG1159 stars with different stellar masses and real- 
istic histories that led them through the thermally puls- 
ing AGB and born again phases (see Dreizler & Heber 
1998; Werner 2001). The lack of a grid of evolutionary 
sequences with different masses in the literature has moti- 
vated us to undertake the present work, which is aimed at 
providing full evolutionary tracks appropriate for PG1159 
stars based on a complete and selfconsistent treatment 
of the progenitor stars, starting from the main sequence. 
Specifically, we follow the evolution of initially 1, 2.2, 2.5, 
3.05, 3.5 and 5.5 Mq stars from the zero-age main se- 
quence (ZAMS) through the thermally pulsing and mass- 
loss phases on the AGB to the PG1159 regime. In partic- 
ular, the evolutionary stages corresponding to the almost 
complete burning of protons following the occurrence of 
the VLTP and the ensuing born-again episode are care- 
fully followed for each sequence. In addition we have fol- 
lowed the evolution of initially 2.2 and 2.7 M model stars 
through the surface hydrogen dilution phase that follows 
the LTP. Attention is paid to the abundance changes dur- 
ing the entire evolution, which are described by means of 
a time-dependent scheme for the simultaneous treatment 
of nuclear evolution and various mixing processes. We also 
develop an evolutionary scenario to address the existence 
of PG1159 stars with unusually high helium abundances 
and nitrogen deficiency, an issue recently raised by Quirion 



et al. (2004) in the context of pulsating PG1159s. We also 
apply this evolutionary scenario to address the existence 
of low-mass O(He) stars. 

The paper is organized as follows. The following sec- 
tion contains the main physical inputs to the models and 
a description of the evolutionary sequences considered. In 
Sect. 3 we present the evolutionary results. There, we elab- 
orate on the main aspects of pre-PG1159 evolution, par- 
ticularly during the born-again phase and the attendant 
chemical changes. In this section we also compare our 
PG1159 models with those existing in the literature. In 
Sect. 4, we discuss the implications of our results. We also 
discuss an evolutionary scenario that could explain the 
existence of PG1159s with high helium content. Finally, 
Sect. 5 is devoted to summarizing our results. 

2. Input physics and evolutionary sequences 

The calculations presented in this work have been done 
with the LPCODE evolutionary code appropriate for com- 
puting the formation and evolution of WD stars through 
late thermal pulses (Althaus et al. 2005a,b). LPCODE 
uses OPAL radiative opacities (including carbon- and 
oxygen-rich compositions) from the compilation of Iglesias 
& Rogers (1996), complemented, at low temperatures, 
with the molecular opacities from Alexander & Ferguson 
(1994). In this work we adopt solar metallicity. The abun- 
dance evolution of 16 chemical elements is included via 
a time-dependent numerical scheme that simultaneously 
treats nuclear evolution and mixing processes due to 
convection, salt fingers and overshooting (OV). Such a 
treatment is particularly important during the short-lived 
phase of the VLTP and the ensuing born-again episode. 
Computationally, convective OV was treated as a expo- 
nentially decaying diffusive process above and below any 
formally convective region, including the convective core 
(main sequence and central helium burning phase), the 
helium flash convection zone and the external convective 
envelope (see Herwig et al. 1999 for details). 

We present six model star sequences, with initially 1, 
2.2, 2.5, 3.05, 3.5, 5.5 M Q . All of the sequences are followed 
from the ZAMS through the thermally pulsing and mass 
loss phases on the AGB (TP- AGB). After experiencing 
several thermal pulses (see Table ^ f° r a description of 
the evolution of the sequences up to the TP- AGB), the 
progenitors depart from the AGB and evolve towards high 
T e g. Mass loss during the departure from the AGB has 
been arbitrarily set, until log(T G ff) ~ 4 when it was turned 
off, as to obtain a final helium shell flash during the early 
WD cooling phase. 

We mention that the 1 M model experienced a core 
helium flash, which led it to the horizontal branch to burn 
helium in a stable fashion. For this sequence two differ- 
ent AGB evolutions have been considered, with different 
mass loss rates as to obtain different number of thermal 
pulses and, eventually, two different remnant masses of 
0.53 and 0.542 M . After mass loss episodes, the stellar 
masses of the post-AGB models is reduced to 0.53, 0.542, 



Miller Bertolami & Althaus: Full evolutionary models for PG1159 stars. 



3 



Initial Mass 




iCHcB 


^TP-AGB 


#TP 


final mass [in M 


1 M 


9790 


12666-12790 


12804 


5/7 


0.530/0.542 


2.2 M 


879 


904-1158 


1183 


16 


0.565 


2.5 M 


595 


609-818 


832 


13 


0.584 


3.05 M 


374 


380-473 


480 


11 


0.609 


3.5 M 


254.5 


257.5-314 


319 


8 


0.664 


5.5 M 


81 


81.6-92.5 


93.3 


8 


0.870 



Table 1. Times for different moments of the evolution of the PG1159 progenitors (in 10 6 yr since the ZAMS). The 
column below icHeB marks the times at the start and the end of core helium burning (CHeB). Also the number of 
thermal pulses is given (column 5). 
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Fig. 1. HRDs for selected VLTP sequences. Legends indicate initial and final masses of the models. In the upper left 
panel, the dotted line corresponds to the case in which a larger number of thermal pulses was allowed before the 
progenitor departs from the AGB. 



0.565, 0.584, 0.609, 0.664 and 0.87. Except for the 1 M Q 
sequences we find convergence problems shortly after the 
departure from the AGB, probably due to the lack of a 
hydrostatic solution in our models (see Wood & Faulkner 
1986); problems that are more noticeable when OPAL 
opacities are used. To overcome this difficulty we have 
temporarily increased the mixing length theory parame- 
ter a by about a factor of two during the early departure 
from the AGB (log T c g <; 3.6). For the purpose of the 
present work, this artificial procedure does not bear any 
relevance for the further stages of evolution. We also want 
to mention that our most massive model (0.870 M Q ) ex- 



perienced several convergence problems, in the envelope 
and atmosphere integration, due to the very high radi- 
ation pressure that characterizes the outer layers of this 
model. This causes the density of the hydrostatic model 
to drop to almost zero and thus our sequence failed to 
converge. To continue with evolution, we imposed arti- 
ficial boundary conditions during its redward excursion 
after the VLTP. We particularly checked that this modi- 
fication did not alter the location of the model in the log 
g-T c ff diagram by doing the same modification to the 3.5 
M Q model, for which the evolution after the VLTP was 
consistently calculated, and we find that no effect of this 
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Fig. 2. Internal chemical profiles for VLTP models hav- 
ing different stellar masses shortly after the star has at- 
tained high luminosities after the VLTP. These chemical 
profiles will be affected by the envelope mixing and mass 
loss episodes when models return to the domain of the 
giant region. Abundances are given in mass fractions. 



artificial evolution was present when the model reached 
the PG1159 phase. 



3. Evolutionary results 

3.1. Evolution during the born-again phase 

In this work we do not intend to present a detailed study 
of the born again phase leading to PG1159 stars, an as- 
pect which will be deferred to a forthcoming publication. 
Instead we limit our discussion to presenting the main 
characteristics of the evolution of the progenitor stars 
relevant for PG1159s. In Fig. 2] we show the complete 
Hcrtzsprung-Russell diagram (HRD) for four of our cal- 
culated sequences (for a complete discussion of the evolu- 
tion of an initially 2.5 M Q sequence see Miller Bertolami 
et al. 2006). As can be noted, our numerical simulation 
covers all the evolutionary phases from the ZAMS to the 
domain of the early- WD including the stages correspond- 
ing to the helium thermal pulses at the tip of the AGB 
and the born again episode. This last stage is particularly 
fast and evolution proceeds on a scale of years (Herwig 
2001, Miller Bertolami et al. 2006). During this stage 
the star is ballooned back to the giant region both by 



the helium flash and by the violent H-burning that oc- 
curs at the VLTP. Note that the detailed evolution of 
the sequences during this stage depends on the mass of 
the remnant. In particular the expansion caused by H- 
burning is very dependent on the mass of the model. We 
find that the double loop path in the HR, characteris- 
tic of successive H- and He-driven expansions (Lawlor & 
MacDonald 2003), is only present for intermediate mass 
models (0.55M Q ^ M <; 0.61M Q ). This is because for 
model with lower masses the H-driven expansion is strong 
enough to keep the model at low effective temperatures 
until the effect of the He-driven expansion takes over, so 
that no return to the blue after the H-driven expansion 
happens. On the other hand, for massive models we find 
that H-driven expansion is extremely weak and it is almost 
suppressed (it stops at log(T e ff) < 4.75), see Fig. 

Another point worth of comment, without entering a 
detailed description of the born again scenario, is that 
the effective temperature at which the model emerges at 
high luminosities is a function of the stellar mass. As can 
be seen in Fig. ^ we find that the lower the mass of the 
model, the lower the effective temperature at which the 
model becomes bright (from log(T e ff) ~ 4 for the 0.53 M Q 
to 5.1 for the 0.664 M ). 

In FigEl we show chemical abundances (abundances 
throughout this article are given in mass fractions) 
throughout the interior of most of the models presented 
in this work (including the previously calculated 2.5 M©) 
at the moment when the models have attained high lu- 
minosities after the VLTP. Note the dependance of the 
composition profile on the stellar mass. The different en- 
velope compositions are due partly to the fact that each 
sequence underwent a different number of thermal pulses 
(see Table 1) during the AGB phase 1 . The abundance pro- 
files shown in FigEl correspond to evolutionary stages well 
before the occurrence of the envelope mixing by convec- 
tion and mass loss episodes that take place at very low 
temperatures after the VLTP. Consequently the depicted 
surface abundances of these models are not strictly the 
ones expected in PG1159s 2 . As can be seen in Fig|2after 
the VLTP models show a chemical structure that can be 
summarized as follows: a C-0 core, an intershell region 
rich in He, C, O and a surrounding envelope where also 
high amounts of 13 C and 14 N (created at the VLTP) are 
found. We elect to show the chemical stratification prior to 
the envelope mixing and mass loss episodes because both 
the deepness of the envelope convection after the VLTP 
and the mass loss rate are highly uncertain and will af- 
fect the final surface composition. For example all of our 



1 We find that the intershell helium abundance reaches a lo- 
cal minimum value after the third-seventh (depending on stel- 
lar mass) thermal pulse (see also Fig. 11 in Herwig 2000). 

2 Because these abundances should be the consequence of 
both the envelope dilution and mass loss. We note that dilu- 
tion of the envelope material with the intershell layers leads 
to lower 14 N surface abundances, between 0.005 and 0.025 by 
mass which is similar to the observed ones (Miksa et al. 2002). 
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Fig. 3. Location of analyzed PG1159 stars ( values taken from Werner & Herwig 2006). Empty symbols are used 
for the objects for which no information about their variability is available (N.D.V.) . Solid lines correspond to our 
evolutionary sequences 0.87, 0.664, 0.609, 0.584, 0.564, 0.542 and 0.53 M©. The broken line corresponds to the (helium 
enriched) 0.512M Q post early AGB model discussed in Section 4.2. 



models remain at log(T eff ) <; 3.8 3 for about 700 - 800 yr. 
Mass loss rates similar to those observed in the Sakurai's 
object (1 x lO _5 M0/yr, and even higher, see Hajduk et 
al. 2005) would cause 7 x 10~ 3 M Q to be stripped off from 
the star 4 . In all the cases the expected surface abundances 
should be somewhat between those displayed by models 
in the envelope and at the intershell region (Fig. EJ) . The 
helium abundance is particularly noteworthy. Indeed note 
that our post-VLTP models predict that the final surface 
abundance of helium will be in the range 0.3-0.55 by mass, 
which is in agreement with the range of observed helium 
abundance in most PG1159s (see Werner & Herwig 2006). 



3 Except for the 0.870 Mq for which we cannot obtain reli- 
able ages. 

4 This strong mass loss at low effective temperatures may 
help to explain the bizarre surface abundance of H 1504+65, 
which is difficult to explain by only considering mass loss dur- 
ing the standard post-AGB evolution (Schonberner & Blocker 
1992). 



However these abundances should be taken with caution 
as they are not only determined by the stellar mass but 
also by the number of pulses considered in the AGB stage 
and the exact value of the adopted OV parameter (/) 
(Herwig 2000). 

3.2. PG1159 stage: Overall location of our models in 
the log g-T eS diagram 

In this section we focus on the location of the models in the 
log g-T c ff diagram. We compare them with the available 
data from observations and from previous stellar models 
for H-dcficient post-AGB stars (Figs. El & El res P ec tfvely) . 
As can be seen in Fig. 01 our models are somewhat hot- 
ter than previous models of He-burning objects with H- 
deficicnt surface abundances. From visual inspection one 
can conclude that the mean shift in spectroscopical masses 
due to this difference in the surface temperatures would be 
0.025 M Q for O'Brien & Kawaler (as presented in Dreizler 
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Fig. 4. Comparison of the models presented in this work 
(thin lines; 0.87, 0.664, 0.609, 0.584, 0.565, 0.542, 0.53 M Q ) 
with previous models of He-burning objects. Dotted lines 
correspond to Wood & Faulkner 1986 type B mass loss 
models (0.87, 0.76, 0.7, 0.6 M Q ) and broken lines indicate 
O'Brien & Kawaler tracks for models with PG1159 like 
compositions (0.63, 0.6, 0.573, 0.55 M Q ). As can be seen 
our models are generally hotter than previous models of 
similar mass. 

& Heber 1998) and 0.04 M for Wood & Faulkner models 
(Wood & Faulkner 1986). On the other hand, we find our 
models to be in good agreement with the 0.604 track pre- 
sented by Werner & Herwig (2006), which is based on sim- 
ilar physical assumptions. This seems to confirm the find- 
ing by these authors that the location of the models in the 
log g-T e ff diagram is very sensitive to the modeling of the 
previous AGB evolution. In particular Werner & Herwig 
(2006) point out that this difference may be related with 
the effect of a smaller effective core growth caused by the 
more efficient third dredge-up events of the new models 
(mainly as consequence of the inclusion of overshooting 
mixing at the He-flash driven convection zone during the 
TP- AGB). If this is true, a lower /-value than adopted in 
the present article (/ = 0.016) at the He-flash driven con- 
vection zone would lead to a less efficient third dredge-up 
and, thus, would shift post-AGB tracks to lower effective 
temperatures 5 . This speculation would be in line with hy- 
drodynamical simulations of the He-shell flash (Herwig et 
al. 2006) that suggest that OV may be reduced to a min- 
imum. However, it should be kept in mind that a quite 

5 After the acceptance of this paper we have not found (in 
preliminary calculations) any important shift for remnants of 
similar mass but different effective core growths during the 
TP- AGB. We intend to explore this effect in a future work 



efficient overshooting at the base of the shell flash driven 
convective zone (during the TP- AGB) is needed in order 
to obtain oxygen surface abundances like those displayed 
by PG1159 stars (Herwig et al. 1999). 

Also regarding the location of post-AGB tracks in the 
log g-Teff diagram, previous works (Werner 2001) have 
noted that it was possible that spectroscopical masses 
which are based on models with normal surface abun- 
dances may suffer from systematic errors. As it is shown in 
Fig. [5] H-burning models are indeed different from those 
which are He-burning of similar mass. However the ef- 
fect is not the same for all masses. We find that differ- 
ences are larger for models of lower mass. In fact our 
H-burning model of 0.53 M Q is located at much lower 
temperatures than the He-burning model of similar mass 
(and thus closer to the location of the low mass models of 
Blocker 1995). Another difference is that in the high-g side 
of the knee, in the log g-T e ff diagram, H-burning models 
are systematically less compact. 

4. Discussion 

4.1. Are PG1159 stars post VLTP or LTP stars? 

It is worth mentioning that all of our VLTP-models dis- 
play high amounts of 14 N at their surface (and also a 
high 13 C/ 12 C ratio see Fig- El ■ This is similar to the re- 
sult of Iben & MacDonald (1995), but at variance with 
later results for the VLTP by Herwig (2001) . Since many 
PG1159s do not exhibit 14 N (Dreizler & Heber 1998), its 
was interesting to analyze other possible channels for the 
formation of PG 1159 stars which do not yield high 14 N 
surface abundances. In particular it was interesting to look 
for changes in the location of the models in the log g- 
T e ff diagram if PG1159 models were the offspring of a 
LTP, in which H is not burnt but diluted after the helium 
flash and thus no enhanced 14 N abundance is expected. 
Consequently we carried out two additional calculations 
of sequences that went through a LTP (of 0.589 and 0.562 
M Q ). These are shown in Fig. where it can be seen their 
location in the log g-T e ff diagram compared with VLTP 
sequences of similar mass. As can be noted there is only 
a slight unimportant difference between VLTP and LTP 
tracks that occurs in the low-gravity region before reach- 
ing the knee 6 . This may be connected to the fact that, 
there, the H-burning shell is still active in LTP sequences. 

In both LTP-sequences, H was diluted to about a 4% in 
mass 7 , which is close to the detection limit (Werner 2001). 
Our LTP-models do not display enhanced 14 N abundance 
at their surface. Our results reinforce the idea of Werner 

6 However while our models are calculated adopting so- 
lar values for the heavy metal abundances, it is known that 
post- VLTP stars display a non-solar heavy metal distribution 
(Asplund et al. 1999 & Miksa et al. 2002). For example, we 
note that tracks move 0.03 dex to the right (red) if opacities 
are increased by ~ 40%. 

7 He/C/O abundances are 31/40/20 and 33/36/21 for the 
0.589 and 0.562 Momodels, respectively. 
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et al. (1998), that the different abundances of 14 N in the 
PG1159 sample seems to be an indication that PG1159 
progenitors are a mixture of VLTP and LTP-AFTP cases. 
As was mentioned by Werner (2001) the existence of LTP 
progenitors in the PG1159 sample may help to solve the 
inconsistency presented by the expansion ages of plane- 
tary nebulae. In this connection we find it interesting to 
discuss the presence of planetary nebulae regarding the lo- 
cation of their central stars in the log g-Tcff diagram and 
their corresponding evolutionary ages. In Fig. 0we draw 
isochrones based on the evolutionary ages of our models. 
We note that ages are measured since the moment of the 
VLTP and not since the first departure from the AGB 8 . 

In order to compare theory with observations, and to 
distinguish between stars coming from VLTP, LTP and 
AFTP scenarios, it seems necessary to see how observable 
parameters like, for example, the presence of planetary 
nebulae, surface abundances and pulsational instability 
correlate with each other. In this context it is interesting 
to see how the presence of high 14 N surface abundances 
correlates with the existence of planetary nebulae around 
PG1159 stars. In what follows we will discuss the impli- 
cations of our evolutionary sequences in the light of the 
correlation between the presence of 14 N and pulsations 
found by Dreizler & Heber (1998). In fact, we will assume 
that this correlation is valid for all PG1159 non-pulsators 
within the instability strip (however see Werner & Herwig 
2006 for a redetermination of nitrogen, 0.001 by mass, in 
the pulsating PG1159-035), even though this correlation 
is based on a very small sample (9 stars). Although this is 
very speculative, we think that inferences may be interest- 
ing enough. Also note that, although 14 N has been demon- 
strated by Quirion et al. (2004) to be unimportant in driv- 
ing pulsations, N is a witness of previous evolution. Stars 
with and without 14 N are expected to have different H 
and heavy elements abundances and consequently differ- 
ent opacities (see Iglesias et al. 1995). As the first travel 
to the high temperature region of the HR (after depart- 
ing from the AGB) takes about 10000-20000 yr (VLTP, 
Blocker 1995) for low mass stars, we find it interesting 
to note that the presence of planetary nebulae around 
some PG1159s close to the isochrones of 25000 and 50000 
yr clearly implies that these objects (particularly NGC 
650, VV47 and PG1520+525) cannot be the offspring of a 
VLTP (because planetary nebulae should have faded dur- 
ing this time interval, 30000-75000 yr). This speculation 
seems to be corroborated by the fact that one of these 
objects (PG1520+525) is known to be 14 N-deficient, and 
also by the fact that the rest of them are non-pulsating 
objects which are close to, or inside, the instability strip 
(and thus according to Dreizler & Heber 1998 must be 
14 N-deficient). It is also worth mentioning that some stars 
(between the 25000 and 100000 yr isochrone) with de- 
tected nitrogen (the ones that should be coming from a 



8 Also note that, as no mass loss has been considered in these 
sequences after the VLTP, these values should be taken as up- 
per limits. 




Fig. 5. Comparison of the location of H and He burn- 
ing post AGB models. Note that at lower g-values (dur- 
ing the WC stage) H-burner are located at slightly larger 
gravities, while at higher gravities (beyond log g ~ 6.5) 
H-burners have lower gravities than He burning objects. 
Also note that the low mass H-burning object lies at lower 
temperatures than the He-burning object of similar mass 
(0.53 M Q ). 

VLTP episode, for instance PG1144+005) do not show an 
associated planetary nebulae, whilst NGC 650 and VV45 
(both non-pulsating objects) do show planetary nebulae. 
This correlation between the presence of 14 N and the ab- 
sence of the planetary nebulae, seems to be another in- 
dicator that 14 N-rich objects should come from a VLTP 
while 14 N-deficient should come from a LTP. 

This argument is not affected by the possibility that 
the pulsating PG1159 at relatively low g (~ 6) (which do 
display planetary nebulae) may be 14 N rich according to 
the correlation found by Dreizler & Heber (1998). In fact, 
because all these objects are more massive than ^0.6 M© 
according to our models, then these objects quickly reach 
the WD cooling track where they experience the VLTP 
(~ 5000 yr, see Blocker 1995). Then, as these objects are 
still close to the 2000 yr isochrone, no more than ~ 7000 
yr have past since their first departure from the AGB and 
thus they can still retain their planetary nebulae, even if 
they have suffered from a VLTP. 

Even beyond the speculation about the correlation be- 
tween 14 N and pulsations, it is interesting to note that 
PG1144+005 and PG1159-035, which do not show asso- 
ciated planetary nebulae (while in the log g-T c g diagram 
they are surrounded by stars that do show associated plan- 
etary nebulae) are known to display 14 N at their surfaces. 
On the other hand, PG1520+525 which is located closer 
to the same isochrone than PG1 159-035 and that is known 



8 



Miller Bertolami & Althaus: Full evolutionary models for PG1159 stars. 




Fig. 6. LTP sequences calculated for this work (0.589 M , 
0.564 M Q , broken lines) are shown agains VLTP sequences 
(thin lines). Also, the helium enhanced post early AGB 
sequence with 0.512 Mq (which is also a post-LTP model) 
is shown (coolest sequence, also in broken lines). Note that 
LTP and VLTP models of similar mass follow almost the 
same track. 

to be 14 N deficient (N< 1 x 10" 4 , Miksa et al. 2002), does 
show a planetary nebulae. 

In this section we mentioned that the existence of 
two groups of 14 N-abundance in PG1159 stars could be 
explained by the existence of two different evolutionary 
channels leading to the PG1159 stage (namely VLTP and 
LTP episodes). If so, our models predict 14 N-deficiency 
to be directly related with the presence of H below the 
detection limit. We also mentioned that this claim might 
be supported by the correlation between the presence of 
14 N and the absence of planetary nebulae. However, it is 
clearly at variance with the reported correlation between 
helium and 14 N mentioned by Quirion et al. (2004), as 
there is no simple reason why LTP models only should 
produce helium-enriched PG1159 stars. In the follow- 
ing section we will discuss an evolutionary scenario for 
these helium-enriched stars (namely HS 1517+7403, HS 
0740+6153 and MCT 0130-1937), an explanation which is 
related to the low mass of these stars. 

4.2. Low T e ff-Low g PG1159 stars. The existence of 
helium-enriched PG1159 stars and their possible 
connection with low mass O(He) stars 

In the previous section we assume the fact that pulsat- 
ing PG1159s are rich in 14 N; however as was noted by 
Dreizler & Heber (1998) and also by Quirion et al. (2004) 




5.5 5.4 5.3 5.2 5.1 5 4.9 4.8 
Log(T ff ) 



Fig. 7. Isochrones for post- VLTP sequences. Times are 
measured since the last helium flash. Empty circles (black 
dots) indicate PG1159 stars with (without) planetary 
nebulae. Note the existence of some objects beyond the 
2.5 x 10 4 yr isochrone that show planetary nebulae. These 
objects could not be the result of a VLTP episode as their 
surrounding nebulae should have faded. Indeed, in none 
of these objects significant amounts of 14 N (a signature 
of previous proton burning during the VLTP) have been 
found while many of them are known to be N-deficient 
(Dreizler & Heber 1998). Also note the fast evolution of 
the most massive model before it reaches log g^ 8, which 
easily explains the absence of detected PG1159 stars with 
log T cff >5.4 

many of the non-pulsator stars in the sample have an 
unusually high helium abundance. This is important in 
view of the results found by Quirion et al. (2004) that 
this helium enrichment would be the cause for the exis- 
tence of non-pulsators within the instability strip 9 . If non- 
pulsators are indeed characterized by helium-enhanced, 
N 14 -deficient surface abundances they cannot be the re- 
sult of the standard LTP scenario experienced by post 
TP-AGB stars. Also, it is clear from visual inspection of 
Fig. |3 that many of the H-deficient stars lie too much to 
the red of the post- VLTP models presented in this work. 
Although comparison with observations has to be taken 
with care due to high uncertainties in the determination 
of g (typical errors in log g are about ~ 0.5 dex), we feel 
that this discrepancy between models and real stars de- 

9 Note, however, that many non-pulsators in the observa- 
tional instability strip are found to have standard helium 
abundances (Longmore 3, Abell 21, VV47, PG1424+535 and 
PG1151-029). On the other hand, NGC 246 which displays a 
high helium abundance (~ 60%, Werner & Herwig 2006) is a 
pulsating PG1159. 



Miller Bertolami & Althaus: Full evolutionary models for PG1159 stars. 



9 




MCT 0130-1937 

o r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 j 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

5.2 5.1 5 4.9 4.8 4.7 4.6 4.5 4.4 

Log(T eff ) 

Fig. 8. Evolution in the log g-T c ff diagram for the 0.512 
M Q remnant. Surface abundances [H/He/ 12 C/ 14 N/ 12 0] 
diplayed by the model are given at certain points (black 
dots) for comparison with the abundances displayed by 
the low mass O(He) stars and the He-enriched PG1159 
stars (triangles and square symbols, respectively). 

serves some discussion. Firstly, it is interesting to note 
that these objects are not probably VLTP descendants, 
as three of them (PG 1424+535, HS0704+6153 and HS 
1517+7403) are known to be 14 N-deficient while no 14 N- 
enhanced PG1159 is known in that region of the log g- 
T e g diagram. Also as it is shown in Fig. [7| some of these 
objects still retain a planetary nebula, a fact that is dif- 
ficult to understand within the VLTP scenario (Werner 
2001). TP-AGB progenitor stars that went through a LTP 
is a possibility to explain the existence of these PG1159s. 
However, as it is shown in Fig. H3 LTP and VLTP mod- 
els of similar masses seem to be located at similar places 
in the log g-T e s diagram. Thus if these stars are descen- 
dents of LTP episodes they should be objects with very 
low stellar masses 10 . 

We have explored an evolutionary scenario that could 
explain the existence of the above mentioned PG1159s. In 
addition we have found that this scenario may offer an al- 
ternative explanation for the formation of low mass O(He) 
stars. This subclass of H-deficient stars is characterized 
by the fact that their members display almost pure He- 
envelopes with only small traces of additional elements. 
They have been proposed to be the offspring of a merging 
event between two white dwarfs (Rauch et al. 2004, see 
also Saio & Jeffery 2002 for the simulation of the merging 
event) and the direct descendants of RCrB stars. However 

10 However it should be noted that stars with very low mass 
may need an unacceptable long time to reach the PG1159 stage 
(< 12.8 Gyr). 



at first glance it would seem that the final masses of the 
products of a merging event (Saio & Jeffery 2002) should 
be higher than the masses derived for the two low mass 
O(He) stars (roughly ~ 0.55M Q ). In what follows we will 
show that the two O(He) stars with low mass (namely 
K 1-27 and HS2209+8229) may be the direct progenitors 
of the He-enriched PG1159 stars. In particular it is worth 
noting that Kl-27, for which detailed abundances have 
been derived, displays a surface composition typical of the 
He buffer of AGB stars, this is almost pure He (~ 98%) 
with traces of 14 N (~ 1.7%) (Werner & Herwig 2006). Also 
both of these stars are supposed to have still an apprecia- 
ble mass loss rate of about ~ 10 _9 MQ/yr, as derived from 
radiation driven wind theory (see for example Table 4 of 
Rauch et al. 1998). 

In order to analyze an evolutionary channel that could 
give rise to 14 N-deficient, He-enhanced PG1159 stars, we 
have calculated one additional sequence for an initially 1 
Mq star, but changing mass loss rates for it to suffer its 
first thermal pulse while it was departing from the AGB 
(so its first thermal pulse is a LTP). This situation is more 
likely to occur in low-mass star progenitors (see Blocker 
1995). The resulting 0.512 M Q remnant is shown in Fig. 
13 and El (post early AGB sequence) . As can be seen this 
model can account for the location of many of the "red- 
der" PG1159 stars. However due to the low intensity of the 
first thermal pulses in low-mass stars no strong dredge-up 
is found and thus the sequence displays a normal surface 
abundance composition after the LTP . We analyzed the 
possibility that the model can turn into a PG1159 due to 
mass loss. This is an important point as the model remains 
in the red giant region (where strong mass loss episodes 
are suppossed to occur) for an unusually long period of 
time after the LTP. To this end we follow the evolution 
using a mass loss rate of 10 -6 M Q /yr if log(T c ff) < 3.8 
and 10" 8 M /yr if log(T cff ) > 3.8. We find that such 
mass loss rates after the LTP can indeed turn the star 
into a PG1159. The fact that the model remains in the 
red giant region for a long period of time is due to the 
several He-subflashes it experiences after the first ther- 
mal pulse (this sub-flashes are typical of the first thermal 
pulses of low-mass stars). It should be noted that, as only 
one thermal pulse has happened, the final He/C/O abun- 
dances of this model (73/21/3) are much different from 
typical PG1159 models. Also as no vigorous H burning 
happened, the star is also 14 N-deficient. In this connec- 
tion it is interesting to note that the three He-enriched 
(and 14 N-deficient) objects in the non-pulsating sample of 
Dreizler & Heber (1998) 11 are located close to this track 
in the log g-T c g diagram and show a very similar com- 
position. Indeed we think that the mentioned correlation 
between He and 14 N may only reflect that the three non- 
pulsators with higher He-enriched surface abundances are 



11 These are HS 0704+6153, HS 1517+7403 and MCT 
0130-1937 (He/C/0= 69/21/9, 84/13/2, 74/22/3 respectively, 
Miksa et al. 2002). 
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Fig. 9. Location of the sequences in the g-T e g-M space. 
Black squares show the derived location for real PG1159 
stars by adopting the T e ff-g values presented by Werner & 
Herwig (2006). 

the offspring of low mass stars that never reached the TP- 
AGB and suffered their first thermal pulse as a LTP. 

In Fig. [H] we show the location of the 0.512 M Q se- 
quence in the log g-T Q g diagram. Additionally the surface 
abundances [H/He/ 12 C/ 14 N/ 12 0] displayed by the model 
are given for certain points during its evolution. The first 
point (at log T e ffPs 4.6 denotes the moment in which sur- 
face abundances start to change as consequence of mass 
loss. At this point layers in which partial H-burning has 
happened are shown by the model. Abundances continue 
changing until the small He-buffer (similar to that of case 
3 in Iben 1995) appears at the surface. At this moment, 
and until the He buffer is eroded by mass loss, the sur- 
face abundance of the models is of 98 % of He and 1.2 
% of 14 N. These abundances are (surprisingly) similar to 
those exhibited by Kl-27, which is located at the same 
region (see FigEJ) like the model in the log g-T c g diagram. 
Finally, when the He-buffer is removed by additional mass 
loss the stars diplays surface abundances typical of the He- 
enriched PG1159. In fact this happens while the model is 
passing through the region (in the log g-T e g diagram) in 
which He-enriched PG1159 are located. Although the ex- 
act locus in the log g-T c ff diagram in which the star will 
show O(He)- or He-enriched- surface abundances depends 
on the exact value of the mass loss rate (which has been 
set somewhat artificially in this work), this numerical ex- 
periment shows that it is possible (with typical mass loss 
rates) that the low mass O(He) stars and the He-enriched- 
PG1159 form an evolutionary sequence. 

5. Derivation of new spectroscopical masses 

By adopting the new tracks we have derived the theoret- 
ical locus of the models in the g-T e ff-M space. We have 
done this by adopting a linear interpolation (i.e. by fit- 
ting a plane) between selected points along the theoret- 



Star 


Previous 


New 


AM 




masses 


masses 




H 1504+65 


0.89 


0.83 


-0.06 


RX J0122.9-7521 


0.72 


0.65 


-0.07 


RX J2117.1+3412 


0.70 


0.72 


0.02 


HE 1429-1209 


0.67 


0.66 


-0.01 


PG 1520+525 


0.67 


0.61 


-0.06 


PG 1144+005 


0.60 


0.55 


-0.05 


Jn 1 


0.60 


0.55 


-0.05 


NGC 246 


0.72 


0.75 


0.03 


PG 1159-035 


0.60 


0.54 


-0.06 


NGC 650 


0.60 


0.54 


-0.06 


Abell 21=Ym29 


0.58 


0.53 


-0.05 


K 1-16 


0.58 


0.54 


-0.04 


Longmore 3 


0.59 


0.54 


-0.05 


PG 1151-029 


0.60 


0.58 


-0.02 


VV 47 


0.59 


0.53 


-0.06 


HS 2324+3944 


0.58 


0.53 


-0.05 


Longmore 4 


0.65 


0.63 


-0.02 


SDSS J001651.42 


0.65 


0.63 


-0.02 


SDSS J102327.41 


0.65 


0.59 


-0.06 


PG 1424+535 


0.57 


0.51 


-0.06 


HS 1517+7403 


0.57 


0.51 


-0.06 


Abell 43 


0.59 


0.53 


-0.06 


NGC 7094 


0.59 


0.53 


-0.06 


SDSS J075540.94 


0.62 


0.58 


-0.04 


SDSS J144734.12 


0.62 


0.58 


-0.04 


SDSS J134341.88 


0.62 


0.58 


-0.04 


SDSS J093546.53 


0.62 


0.58 


-0.04 


SDSS J121523.09 


0.62 


0.58 


-0.04 


HS 0444+0453 


0.59 


0.55 


-0.04 


IW 1 


0.56 


0.50* 


-0.06 


Sh 2-68 


0.55 


0.50* 


-0.06 


PG 2131+066 


0.58 


0.55 


-0.03 


MCT 0130-1937 


0.60 


0.54 


-0.06 


SDSS J034917.41 


0.60 


0.54 


-0.06 


PG 1707+427 


0.59 


0.53 


-0.06 


PG 0122+200 


0.58 


0.53 


-0.05 


HS 0704+6153 


0.51 


0.47* 


-0.04 



Table 2. Masses derived from the new tracks presented in 
this work by adopting the values of T c g and g from Werner 
& Herwig (2006). Previous masses, which are based on 
old post-AGB tracks, have been taken from Werner & 
Herwig 2006. Asterisks denote masses that have been de- 
rived by linear extrapolation and that are consequently 
more uncertain than the others. The new mean mass for 
the PG1159 sample is O.573M , 0.044M o lower than pre- 
vious mean mass. All masses are in Mq. 



ical tracks 12 . Once the theoretical surface on the g-T e g- 
M space was derived we have calculated spectroscopical 
masses by adopting the T e g-g values presented by Werner 
& Herwig (2006). In Fig. El we show both the theoretical 
surface and the derived parameters for real PG1159 stars. 
The resulting spectroscopical masses are shown in Table 

m 



Previous to this we tried a global least square polinomical 
fitting but it turned to be unappropriate. 



Miller Bertolami & Althaus: Full evolutionary models for PG1159 stars. 



11 



6. Conclusions 

In this paper we have presented full evolutionary calcu- 
lations appropriate to post-AGB PG1159 stars. In our 
calculations, the complete evolutionary stages of PG1159 
progenitors having a wide range of initial stellar masses 
have been considered, particularly the evolution through 
the born again stage and late thermal pulse (LTP) that 
occur following the departure from the AGB. The loca- 
tion of our PG1159 tracks in the log g-T c ff diagram and 
their comparison with previous calculations as well as the 
resulting surface composition have been discussed at some 
length. In particular, our models are generally hotter than 
previous models of similar masses which do not display an 
efficient third dredge up during the TP-AGB, but similar 
to recent models that include OV (and in which efficient 
third dredge up is found during the TP-AGB). This rein- 
forces the idea that theoretical log g-T e ff diagram is sen- 
sitive to the previous evolution on the AGB. 

Our results reinforces also the idea that the differ- 
ent abundances of 14 N observed in the surface of those 
PG1159 stars with undetected hydrogen is an indication 
that the progenitors of these stars would have evolved 
through born again episodes, where most of the hydrogen 
content of the remnant is burnt, or late thermal pulses 
where hydrogen is not burnt but instead diluted to very 
low surface abundances. 

We have also derived isochrones from our PG1159 
models aimed at discussing the presence of planetary neb- 
ulae in connection with the location of their central stars 
in the log g-T e ff diagram and their corresponding evolu- 
tionary ages. In this regard, we infer that the progeni- 
tors of some PG1159 stars (particularly NGC 650, VV47, 
IW, and PG1520+525) might not have experienced a born 
again episode, in agreement with the correlation between 
the presence of 14 N and pulsating PG1159s (Dreizler & 
Heber 1998). We also discuss the correlation between the 
presence of planetary nebulae and the 14 N abundance as 
another indicator that 14 N-rich objects should come from 
a born again episode while 14 N-dcficient should come from 
a LTP. 

Also, based on the new models presented in this ar- 
ticle, we derived new values for spectroscopical masses. 
We find that new masses are systematically lower and 
consequently at variance with those derived from aster- 
oseismology. Although typical uncertainties in spectro- 
scopical masses, and also uncertainties in asteroseismolog- 
ical masses (as they are model- and method- dependent, 
Quirion et al. 2004, Corsico & Althaus 2006), prevent us 
from making concluding remarks, we think that this dif- 
ference may be indicating that the microphysics or the 
previous evolution (as for example the number of ther- 
mal pulses on the AGB) may be somewhat different than 
adopted in this work. Both improved determinations in g 
and T ff for PG1159 stars and the use of realistic PG1159 
models in asteroseismological analysis are in order to solve 
this problem. 



Finally, we discuss the possibility that the He-enriched 
PG1159 stars can be represented by very low mass ob- 
jects, the progenitors of which have abandoned the AGB 
shortly before the occurrence of the first thermal pulse. 
These progenitors experience their first thermal pulse as 
a LTP. We find that moderate mass loss rates operating 
during the evolution through red giant region after the 
LTP could turn these objects into H-deficient ones, but 
with high surface helium abundances and deficiency in ni- 
trogen. In particular, we find that the surface abundances 
of HS 0704+6153, HS 1517+7403 and MCT 0130-1937, 
three helium enriched objects in the sample of Dreizler 
& Heber (1998), can be reproduced within this scenario. 
Also our numerical experiment shows that this scenario 
provides a possible explanation for the existence of O(He) 
stars of low mass, linking them as the immediate ancestors 
of the He-enriched PG1159. 

Finally, tabulations of T e ff, L+ and g for the tracks 
presented in this work are freely available at our 
URL;http : //www. f caglp .unlp . edu. ar/evolgroup 
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